Mapping Image Potential States on Graphene Quantum Dots 
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Free electron hke image potential states are observed in scanning tunneling spectroscopy on 
graphene quantum dots on Ir(lll) acting as potential wells. The spectrum strongly depends on 
the size of the nanostructure as well as on the spatial position on top, indicating lateral confine- 
ment. Further analysis of the substructure of the first state by spatial mapping of constant energy 
local density of states reveals characteristic patterns of confined states. The most pronounced state 
in the spectra is not the ground state, but an excited state with a favorable combination of the local 
density of states and parallel momentum transfer in the tunneling process. Chemical gating tunes 
the confining potential by changing the local workfunction. Our experimental determination of this 
workfunction allows to deduce the associated shift of the Dirac point. 



PACS numbers: 73.22.-f, 73.63.Hs, 73.22.Pr, 73.20.-r 

Confinement of electrons in nanostructures leads to 
quantum size effects as a size-dependent electronic struc- 
ture and the appearance of atom-like states (character- 
ized by a set of quantum numbers). In recent years, first 
experiments regarding the confinement of image potential 
states (IPSs) using the spatial resolution of the scanning 
tunneling microscope (STM) have been performed [1 6j, 
transcending the pioneering studies based on two pho- 
ton photoemission (2PPE) [71 [8]. IPSs are unoccupied 
states in an attractive image charge Coulomb potential 
between the Fermi and the vacuum level. Perpendicu- 
lar to the surface they feature a hydrogen-like spectrum 
^9j, whereas parallel to it a two dimensional electron gas 
(2DEG) forms. In STM, IPSs appear as peaks in the lo- 
cal density of states (LDOS) measured while retracting 
the tip from the surface. As they are Stark-shifted due 
to the electric field between tip and sample [lOj they are 
often referred to as field emission resonances (FERs). 

In order to observe and quantify confinement effects 
for IPSs, the corresponding nanostructures must fulfill 
four conditions: (i) a local workfunction ^ smaller than 
that of the surrounding area [2j, a large difference pro- 
vides strong confinement, (ii) a well-defined shape, (iii) 
an established preparation that allows to adjust the size 
in a wide range, and (iv) stability under the high bias 
voltage used in the STM measurements. Whereas pre- 
vious work provides fascinating first insights into quan- 
tum size effects for IPSs, no study yet matches all four 
conditions: In [3j a first hint at a size dependence of 
the energies of IPSs confined to Co islands on Au(lll) 
is visible, however here the size variation was less than 
an order of magnitude. Atom-like patterns have been 
observed above stacking- fault tetrahedra on Ag(lll) [5j, 
still the difference in ^ is so small that the resulting weak 
confinement only acts on the IPSs lowest in energy. The 
system NaCl on metal is promising as it shows a large 
difference in However, up to now there is no estab- 
lished method to tune the size of islands with a well- 



defined shape over a wide range [TP,T2]. In consequence, 
in these experiments electron confinement has not been 
observed yet. Strong confinement is found for islands of 
alkali metals on Cu(lOO) In this case, the atomic 
structure of the islands is unclear, the size cannot be 
varied, and the clusters are not entirely stable during the 
measurement. A feature of special interest is the cou- 
pling between the IPSs on neighboring nanostructures to 
molecule-like states [3| [6l [T3]. 

In this study, we demonstrate that confinement of IPSs 
can be observed in nanostructured epitaxial graphene 
(gr). IPSs on graphene are of special interest as they 
reveal specific properties of this material: As a conse- 
quence of its 2D-character a splitting of the first IPS has 
been predicted [14J and observed for gr/SiC [15j. Split- 
ting of IPSs has also been observed for gr/Ru(0001) due 
to the strong corrugation of the carbon sheet which al- 
lows a large probability density also between graphene 
and the substrate p!6[ ITTJ . In the system under investiga- 
tion here, 2PPE could demonstrate the parabolic bands 
of IPSs in the large band gap of the Ir substrate [18]. 
However, no effect of confinement could be determined 
due to the averaging nature of this method. On funda- 
mental grounds, the IPSs on graphene share a common 
origin with specific electronic states of related materials 
[14J as the interlayer state of graphite, or superatomic 
states of fullerenes [13j. 

For epitaxial graphene, the confining potential can be 
modified: The intercalation of electron acceptors (as, e.g. 
O [19j) between the carbon sheet and its substrate leads 
to a depletion of the charge density in graphene's tt- 
system, which in turn shifts the Dirac point to higher 
electron energies [20j, and vice versa for electron donors. 
Note that this chemical gating of graphene allows signif- 
icantly larger shifts of £^d than electrostatic gating of- 
ten applied in a field effect transistor geometry. Under 
the assumption that the band structure of graphene is 
rigidly pinned to the vacuum level [21j, this shift is di- 
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rectly related to the change in the workfunction of the 
intercalated graphene ^gr/x with respect to the pristine 
case ^gri $gr/x — ^gr = ^D- In consequence, the local 
doping level can be determined via a measurement of 
the local workfunction. 

Nanostructured graphene is found in the system 
graphene quantum dots (GQD) on Ir(lll) studied in 
this work, which fulfill all conditions outlined above: 
(i) a large workfunction difference [A<l> = <l>gr — $ir = 
(4.65 ± 0.10) eV - (5.79 ± 0.10) eV = (-1.1 ± 0.1) eV 
[18]], (ii) a well defined polygonal shape which can be 
determined with atomic precision [22], (iii) a size tun- 
able from less than 10 nm^ to electronically equivalent 
to infinite ^23j, and (iv) high stability also for high bias 
voltages due to the strong carbon-carbon bonds as well 
as the good conductivity of graphene. Fullerenes [13j and 
carbon nanotubes [3j can be viewed as extreme cases of 
confinement. However, for such systems a tuning of the 
size over orders of magnitude is impossible and the curved 
geometry gives a new character to the now strongly hy- 
bridized states. 

Our Ir(lll) substrates are cleaned in ultra high vac- 
uum by cycles of 1.5 keV Ar+ ion bombardment at room 
temperature, oxygen firing at 1120 K and annealing at 
1470 K. Graphene is prepared by decomposition of ethy- 
lene in a temperature programmed growth process at 
1270 K, resulting in a coverage of about 22 % of a mono- 
layer [23j. The STM tip is virtually grounded and the 
sample is put to the bias voltage U leading to a tunneling 
current /. Energies are connected to U by E — E^ = e-U. 
IPSs are investigated by measuring the differential con- 
ductivity in form of both dI/dU{E — E-p) point spec- 
tra and constant energy maps in constant current mode 
(stabilization values t/gtab, ^stab) with active feedback 
loop by using the lock-in technique (/ = 1.317 kHz, 
^mod = 14 mV). Each point spectrum represents the av- 
erage of at least two single measurements. Active feed- 
back is of vital importance as moving the tip in z direc- 
tion during point spectroscopy compensates for the Stark 
shift within one measurement [24J. The tip traces z{V) 
were recorded parallel to the spectra. All data is taken at 
5.3 K in a background pressure lower than 1 x 10~^^ mbar. 
Data analysis is performed using WSxM [25J. Fig.[l](a) 
shows a set of GQDs with sizes ranging from less than 
10 nm^ up to several hundred nm^. A majority of hexag- 
onally shaped islands coexists with irregular ones result- 
ing from coalescence [23j. The d//d[/-spectra taken with 
the same microscopical tip [26j at the center of the flakes 
labeled in Fig.[l](a) show very pronounced features, see 
Fig.[l](b) for the smallest and the largest hexagonal flake. 
We attribute the discrete peaks to the energies E'^^^ of 
the sequence of IPSs with order n = 1, 2, .... These ener- 
gies are down-shifted with respect to pristine Ir(lll) due 
to the significantly smaller ^. However, the spectra on 
the differently sized flakes are also shifted with respect 
to each other by as much as AE ^0.67 eV between the 
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FIG. 1. (color online), (a) Set of graphene quantum dots on 
Ir(lll) [some of them labeled (A)-(G) with increasing size]; 
image size 55nmx90nm, [/ = IV, / = 0.2nA. {h)dI/dU 
spectra on flakes (A) and (E); [/stab = 1 V, /stab = 0.2 nA. 



largest and the smallest flake for both E'^'^^ and E'^^\ Un- 
der our experimental conditions, such a large difference 
cannot be due to a variation of the Stark shift [27J. In 
addition to the size-dependence of the peak energy, the 
individual peaks show a clear substructure, especially the 
ones for the smaller island (red line). This is in contrast 
to the smooth peaks observed for extended systems. 

A dependence of the energy of electronic states on the 
size of the system as well as the occurrence of discrete 
eigenstates evidenced by the peak substructure is a clear 
fingerprint of lateral quantum confinement. To allow an 
analytic treatment of our system, we approximate our 
hexagonal GQDs by an infinite cylindrical potential well 
with radius r [28j. The eigenfunctions in polar coordi- 
nates (p, (/)) are then "^rn^i ^ JiiJ^m.ip)^^^^^ i where Ji is 
the spherical Bessel function of the first kind with order / 
and km,i is the momentum of the electron [29j. Due to the 
confinement, ^rn,i must have a node for p = leading to 
the condition km,ir = Zm,i with Zm,i the m-th zero of Ji. 
Hence, the eigenstates can be labeled (m, /). Radial cuts 
through the normalized probability density ^^m,^ of 
the first six eigenstates are shown in Fig.[2](a). Their en- 
ergies are given by E^^j = E^^^ + (n\/2) • zl^^Ji^l ■ m*). 



3 



(a) 
2.0 

1.5 

CO 

§1.0 
0.5 
O.Q 



' 1 
N3,0) 


1 1 1 1 1 1 1 


1 ' 1 ' 1 ' 
/c||''°^=1.56 1/nrTT 

/c[| ''^=2.48 1/nm' 


^\ 




/cj'°^=3.57 1/nm- 
/c|2'^M.54 1/nm" 


"(1,0) \ 


(2,2) 


/cf''^=5.45 1/nm. 
/cf|^°^=5.60 1/nm- 



.2 0.4 



0.6 0.8 
r (nm) 



1.0 1.2 1.4 




0.00 0.02 0.04 0.06 0.08 0.10 0.12 
inverse area 1/Q [Mnnf] 



0.14 



FIG. 2. (color online), (a) Normalized LDOS and k\\ for a 
model island with Q = 7.5 nm^, i.e. r = 1.545 nm (right bor- 
der of the plot), (b) Energy positions for n = 2 and n = 3 
depending on inverse area (black squares) , measured on quan- 
tum dots (A)-(G) in Fig.[l](a) with parameters from Fig.[l](b); 
peak FWHM (vertical bars) and linear fits (black lines); cal- 
culated energies for (m,/), see labels. 



with E^""^ the energy of the state on extended graphene, 
Q the area of the flake and m* the effective electron mass. 

In Fig.[2](b) we show the dependence of E^'^'> and E^^'' 
on for the whole set of flakes shown in Fig.[l](a). The 
bars correspond to the Gaussian full width half maximum 
(FWHM) of the respective peaks. A linear dependence 
can be fitted to the data (thin black lines), as expected for 
^mi' yielding Eq^^ as the intercept. With m* ~ rrie [18j, 

we can calculate E^\. As we have measured the spectra 
always in the center of the flakes [see also Fig.[l](a)], we 
expect the dominant contribution from eigenstates with 
non vanishing amplitude at the center, i.e. / = 0, marked 
by the labeled lines in Fig.[2](b). Obviously the data fits 
best to the (2,0)-state, both for n = 2 and n = 3. This 
is surprising as this is not the ground state of the sys- 
tem. The dominance of (2,0) cannot be explained by 
the LDOS of the different eigenstates at the center of 
the flake, as this quantity increases monotonically with 
m for / = [see Fig.[2](a)]. For a correct interpreta- 
tion, one has to take the details of the tunneling process 



into account: Most of the tunneling electrons have a van- 
ishing momentum parallel to the surface. However, for 
the electrons confined above the GQDs, km,i increases 
with increasing order [see the values in Fig.[2](a)], mak- 
ing them less accessible for STS. In consequence, whereas 
the LDOS in the center increases with the order of the 
state, the contribution to the tunneling density of states 
(TDOS, e.g. [30J) measured in STM decreases. This 
fact can be quantified by TDOS = LDOS • exp(— zq/A) 

[TlJ, with = 2^2me^/h^ ^k^^ i and k^^i = Zm,i/r. 
Therefore, the contribution of a given state (m, /) to the 
signal observed in STM depends on the tip-sample dis- 
tance zo, the radius of the flake r, and the order of the 
state. Evaluating the above formula for the model is- 
land with Q = 7.5 nm^ indeed yields that (2,0) domi- 
nates for 0.36 nm< zq < 1.04 nm. It is reasonable to 
assume that zq in our experiment is within these bound- 
aries. The considerations on the tunneling probability 
and thus preferred states for tunneling may also connect 
to previous publications like [32j, where the observed de- 
crease in peak intensity could also be explained by less 
probable tunneling due to increased k\\ instead of an ad 
hoc assumption of a peak broadening increasing with en- 
ergy. 

The spatial modulation of the LDOS can be resolved in 
two dimensions by dl/dU mapping on a hexagonal flake, 
see Fig.|3l Here, we chose n = 1 in order to avoid the 
experimental difficulty of obtaining a high spatial resolu- 
tion at high bias voltage (the particle-in-a-box behavior 
of IPSs does not depend on n [3\). Again, the maximum 
of the TDOS-curve corresponds to (2,0), see Fig.[3](f). 
Indeed, the maps corresponding to the selected energies 
bear a strong resemblance to the eigenstates as depicted 
in Fig. [2] (a): The state shown in Fig.jsja) shows a broad 
maximum resembling (1,0), (b) is more peaked in the 
center like (2,0), and (c) and (d) have vanishing intensity 
in the center like (2,1). Strictly speaking, however, the 
patterns we observe are not pure states, but a superpo- 
sition of several neighboring states. Note that in Fig. [3] 
(d) and even more for (e) the breaking of the cylindrical 
symmetry by the hexagonal shape of the flake becomes 
evident. Similar patterns as the above have been ob- 
served on GQDs as a result of confinement of low energy 
occupied states [22j. Finally, the inversion of contrast 
in Fig.[3](e) indicates that the LDOS at this energy is 
dominated by Ir IPSs, drawing our attention to the in- 
teraction of the IPS-2DEGs above graphene and Ir(lll). 
The whole series of IPSs is down-shifted on a GQD with 
respect to Ir(lll) due to the lowered ^. Mapping the 
IPSs along lines cutting through a GQD (see [27J) shows 
that this shift for n = 1 across the boundary region is 
rather abrupt, whereas all higher orders show a more 
continuous way of change. This indicates a suppressed 
interaction for n = 1 and interacting IPS-2DEGs of the 
flake and the substrate for n > 1. 
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FIG. 3. (color online). (a)-(e)IPS n = 1 by 2D constant 
energy mapping of dl/dU on a Q = 11 nm^ sized gr island. 
(a)-(e): E - Ef = 4.50 eV, 4.73 eV, 4.82 eV, 4.98 eV, 5.22 eV; 
/ = 0.2nA; image size: all 5.7 x 5.7nm^ (i) dl /dU{E - Ef) 
spectrum of n = 1 with energies of maps shown in (a)-(e) 
(black squares), E^\ indicated by solid (dashed) lines for / = 

(1) oc ?7^(2>3) Fig.[2](b). 



(1^0); E^Q^ as E^^^'^^ in 



In the following, we will exemplify chemical gating of 
graphene for the case of O. The preparation described 
above is extended by exposing the sample to 750 L of 
O2 at 430 K. This process leads to O intercalation for 
all but the smallest flakes [19j. In order to demon- 
strate the effect of intercalation most clearly, we will 
focus on the largest flakes which are representative for 
extended graphene. Fig. |4] (a) shows an STM image 
of 0-intercalated graphene. The superstructure already 
visible in this image can be accurately determined as 
(V3 X V^)R30° usine; the Fourier transform shown in 
Fig. [4] (c), i.e. a coverage of 0.33 ML. In Fig.ji] (b) we 
compare spectra on gr/O/Ir (blue line) and on 0/Ir(lll) 
(black line). The difference in ^ can be determined from 
a plot A^;^^) = £^gr/o/ir - ^o/ir ^^^^^^ n [331, see in- 
set of Fig.|4] (b), as A<l> = (-1.3 ± 0.1) eV, leading to 
^gr/o/ir = (5.1 ±0.1) eV, which has to be compared with 
^gr/ir = (4.7±0.1)eV HaEJ. 

According to [21], one can deduce A£^d = ^D,gr/o/ir — 
^gr/ir = (0.4 ± 0.1) eV from our 
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experiments. This nicely agrees with A^d = 0.3 eV 
implied by a recent photoemission study [5U] (linearly 
interpolated using the experimentally determined value 
of A£^D = 0.6 eV for 0.6 ML intercalated O). In conse- 
quence, our determination of the local workfunction pro- 
vides direct access to the position of the Dirac point £^d, 
which is often hard to determine by other methods. As an 
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FIG. 4. (color online), (a) topography: O intercalated GQD, 
size 7.6 X 7.6 nm^; (b) IPS spectra on gr/0/Ir(lll) (blue) 
and on 0/Ir(lll) (black), A$ indicated by black arrow; in- 
set: AS(») = - 4%; (c) FT-STM of gr/0/Ir(lll), 

squares: F-points of graphene, circles: (v^ x v^)R30°-O in- 
tercalation superstructure. 



example, for gr/Ir(lll) the LDOS determined using STS 
does not show a pronounced dip at [35j. A determi- 
nation of via ^ can be especially useful for mapping 
the doping level in inhomogeneous graphene systems. 

As an outlook, the suitability of graphene as build- 
ing material for the nanostructures as demonstrated here 
opens new possibilities for the investigation of laterally 
confined IPSs, as there are various methods on how 
to prepare rationally designed architectures based on 
graphene: It can be cut by STM lithography [36], allow- 
ing complicated well geometries. A clever choice of the 
hydrocarbon precursor leads to the formation of superpe- 
riodic structures f37] which may give rise to backfolding 
of the IPS bands. One can envision that it is possible to 
move GQDs with an STM tip moving in close proximity 
to the surface, allowing exact studies of the interaction 
of neighboring quantum wells. Such structures would 
resemble a diatomic molecule, where the participating 
quantum wells can even have different energy levels due 
to a variation in size or in doping. 

Summing up our results, the large difference in work- 
function between graphene and Ir(lll), the well defined 
shape of the nanostructures as well as their large size 
variation and high stability enabled us to demonstrate 
confinement effects of IPSs. We have shown that the 
energy spectrum depends on the size of the GQDs and 
evolves into a series of atom-like states, which is dom- 
inated by a state other than the ground state due to 
an interplay of density of states and parallel momentum 
transfer in the tunneling process. Intercalating extended 
graphene with an electron acceptor introduces an addi- 
tional degree of freedom as this allows tuning of ^. The 
determination of the local workfunction allows to deduce 
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the local doping level £^d- Our findings open new possi- 
bilities for the study of quantum size effects as graphene 
is a very flexible building material for nanostructures. 
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